Integrable theory of quantum transport in chaotic cavities 
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The problem of quantum transport in chaotic cavities with broken time-reversal symmetry is 
shown to be completely integrable in the universal limit. This observation is utilised to determine 
the cumulants and the distribution function of conductance for a cavity with ideal leads supporting 
an arbitrary number n of propagating modes. Expressed in terms of solutions to the fifth Painleve 
transcendent and/or the Toda lattice equation, the conductance distribution is further analysed in 
the large-n limit that reveals long exponential tails in the otherwise Gaussian curve. 
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Introduction. — The low temperature electronic con- 
duction through a cavity exhibiting chaotic classical dy- 
namics is governed by quantum phase-coherence effects 
1 , ]2j . In the absence of electron-electron interactions 
3 , 0, HI , the most comprehensive theoretical framework 
by which the phase coherent electron transport can be 
explored is provided by the scattering S-matrix approach 
pioneered by Landauer There exist two different, 
though mutually overlapping, scattering-matrix descrip- 
tions [3] of quantum transport. 

A semiclassical formulation 0] of the S-matrix ap- 
proach is tailor-made to the analysis of energy-averaged 
charge conduction [9( through an individual cavity. Rep- 
resenting quantum transport observables (such as con- 
ductance, shot-noise power, transferred charge etc.) in 
terms of classical trajectories connecting the leads at- 
tached to a cavity, the semiclassical approac h l lOl effi- 
ciently accounts for system-specific features [11| of the 
quantum transport. Besides, it also covers the long-time 
scale universal transport regime 12| emerging in the limit 
[IH td > te, where td is the average electron dwell time 
and te is the Ehrenfest time (the time scale where quan- 
tum effects set in). 

The latter universal regime [3] can alternatively be 
studied within a stochastic a ppr oach 0, HH based on 
a random matrix description [19 ] of electron dynamics 
in a cavity. Modelling a single electron Hamiltonian by 
an M x M random matrix !K of proper symmetry, the 
stochastic approach starts with the Hamiltonian H to t of 
the total system comprised by the cavity and the leads: 



Htot = X) ^fc M «^ + X 
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Here, xpk and \a are the annihilation operators of elec- 
trons in the cavity and in the leads, respectively. In- 
dices k and I enumerate electron states in the cavity: 
1 < k,£ < M, with M — * oo. Index a counts propa- 
gating modes in the left (1 < a < Nl) and the right 



(N L + 1 < a < N) lead. The M x N matrix W describes 
the coupling of electron states with the Fermi energy sf 
in the cavity to those in the leads; N = TVl + is the 
total number of propagating modes (channels). Since in 
Landauer-type theories the transport observables are ex- 
pressed in terms of the N x N scattering matrix [f| 

S(e F ) = 1 - 2i 7 rW t (e F - JC+ i7rWW t ) _1 'W, (2) 

the knowledge of its distribution is central to the stochas- 
tic approach. (Two such observables - the conductance 
G = tr(CiSC2§^) and the shot noise power P = 
tr (Ci§C2S t ) — tr (ei§62§^) 2 measured in proper dimen- 
sionless units [3] - are of most interest. Here, Ci = 
diag(ljv L , 0at r ) and 62 = diag(0Ar L , 1jv h ) are the pro- 
jection matrices.) 

For random matrices !K drawn from rotationally invari- 
ant Gaussian ensembles [Tjj, the distribution of §(ef) is 
described [HI by the Poisson kernel [3, [3, H(| 



P(S) oc [det (l - SSt) det (l - SS 1 ')] 



P/2-1—PN/2 



(3) 



Here, (3 is the Dyson index (17J accommodating sys- 
tem symmetries {(3 = 1, 2, and 4) whilst § is the av- 
erage scattering matrix S = diag(^/l — FJ) V, 
that characterises couplings between the cavity and the 



leads in terms of tunnel probabilities [21| Tj of j-th 
mode in the leads (1 < j < N); the matrix V is 
V G G(N)/G(N L ) x G(N R ) where G stands for orthog- 
onal {(3 = 1), unitary ((3 — 2) or symplectic ((3 = 4) 
group. 

The above description becomes particularly simple for 
chaotic cavities that coupled to the leads through ballis- 
tic point contacts ("ideal" leads, Tj = 1). Indeed, unifor- 
mity of P(§) over G(N) implies that scattering matrices 
S belong [22j to one of the three Dyson circular ensembles 
[13] about which virtually everything is known. Notwith- 
standing this remarkable simplicity, available analytic re- 
sults for statistics of electron transport are quite limited 
0, IIH . In particular, distribution functions of conduc- 
tance and shot noise power, as well as their higher order 
cumulants, are largely unknown for an arbitrary number 



2 



of propagating modes, TVl and ATr, and thus do not catch 
up with existing experimental capabilities 24 1 . 

In this Letter, we combine a stochastic version of the 
S-matrix approach with ideas of integrability 25, 26| to 
show that the problem of universal quantum transport 
in chaotic cavities with broken time-reversal symmetry 
(/3 = 2) is completely integrable. Although our theory 
applies j2?| to a variety of transport observables, the fur- 
ther discussion is purposely restricted to the statistics of 
Landauer conductance. This will help us keep the pre- 
sentation as transparent as possible. 

Conductance distribution. — In order to describe fluc- 
tuations of the conductance G = tr (Ci§C2§^) in an ade- 
quate way, one needs to know its entire distribution func- 
tion. To determine the latter, we define the moment gen- 
crating function 

J n {z) = (exp (~zG)) SeCVE(2n+l/) (4) 

which, in accordance with the above discussion, involves 
averaging over scattering matrices § € CUE(2n + v) 
drawn from the Dyson circular unitary ensemble (l7j . For 
the sake of convenience, we have introduced the notation 
n = min(A r L, N-r) and v — \N\^ — so that the to- 
tal number TVl + ATr. of propagating modes in two leads 
equals 2 71 + v. 

While the averaging in Eq. ^ can explicitly be per- 
formed with the help of the Itzykson-Zuber formula [28j |. 
a high spectral degeneracy of the projection matrices Ci 
and 62 makes this calculation quite tedious. To avoid 
unnecessary technical complications, it is beneficial to 
employ a polar decomposition (l9| of the scattering ma- 
trix. This brings into play a set of n transmission eigen- 
values T = (Ti,-- - ,T n ) € (0,1)" which characterise 
the conductance Q in a particularly simple manner, 

g(t) v; ; r,. ' 

The uniformity of the scattering S-matrix distribution 
gives rise to a nontrivial joint probability de nsity function 
of transmission eigenvalues in the form [29j, 1301 ] 



with 



P n {T) = c- 1 A 2 n (T 



n^r ( 5 ) 

is the Vandcrmonde de- 



Here A n (T) = Uj <k (T k - Tj _ 
tcrminant and c„ is the normalisation constant 
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Let us stress that the description based on Eq. (0 is com- 
pletely equivalent to the original, microscopically moti- 
vated 8 € CUE(2n + v) model. 

Now the moment generating function can elegantly be 
calculated. A close inspection of the integral 



n. n 



(7) 



reveals that it admits the Hankel determinant represen- 
tation [H 



Jr(z) 



("+!)! 




(9) 



In deriving Eqs. (|SJ) and (J9]) we have used the Andreief-de 
Bruijn integration formula [3lj |. 

Equation ([5]), supplemented by the "initial condition" 
$o(z) — 1, has far-reaching consequences. Indeed, by 
virtue of the Darboux theorem 32j , the infinite sequence 
of the moment generating functions (3^, 3^, • ■ ■ ) obeys 
the Toda lattice equation (n > 1) 

J n (z) 3%(z) - (J' n (z)f = var„(G) J n _i(z) ? n+1 (z), 

(10) 

where var n (G) = n(n + l) 1 (c„_ic, i+ i/c^) is nothing but 
the conductance variance 



var„(G) 



n 2 (n + v) 2 



(2n + v) 2 [{2n + v) 2 - 1] ' 



(11) 



Since 3^ n (z) is the Laplace transform of conductance 
probability density f n {g) = (8{g — GYh the Toda lattice 
equation provides an exact solution [331 ] to the problem 
of conductance distribution in chaotic cavities with an 
arbitrary number of channels in the leads. Equations © 
- (fill) represent the first main result of the Letter. 

There exists yet another way to describe the conduc- 
tance distribution. Spotting that the moment genera ting 
function 3r n (z) is essentially a Fredholm determinant 34 1 
associated with a gap formation probability 17] within 
the interval (z, +00) in the spectrum of an auxiliary nxn 
Laguerre unitary ensemble, 



,-A.f 



A*(A), (12) 



one immediately derives 3, 35|: 

ov(*) - n(n + v) 



J„(z) = exp ( J (II- 



(13) 



Here, ay(t) satisfies the Jimbo-Miwa-Okamoto form of 
the Painleve V equation [36| 

(t<r v ) 2 + Wv - to-' v + 2{a' v ) 2 + (2n + u)a' v ] 2 

+ 4^) Vv + n)(*v + n + v) = (14) 

subject to the boundary condition ay(t — > 0) ~ n(n+p). 

To the best of our knowledge, this is the first ever ap- 
pearance of Painleve transcendents in problems of quan- 
tum transport. The representation Eq. (I13p . being the 
second main result of the Letter, opens a way for a non- 
perturbative calculation of conductance cumulants. 



J n (z) = -det [{-d z ) j+k ?!(*)] 



(8) 
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Conductance cumulants. — Our third main result is the bilinear recurrence relation (j > 2) 
l(2n + v) 2 -j 2 ](j + l) K]+1 = 2^(3^+l)(j-£) 2 ( 3 J - (2n + i/)(2j - 1) jk, - l)(j - 2) (15) 



for conductance cumulants {kj}. Taken together with 
the initial conditions provided by the average conduc- 
tance Ki = n(n + v)j (2n + v) and the conductance 
variance K2 = n\/[(2n + v) 2 — 1], this recurrence ef- 
ficiently generates (previously unavailable) conductance 
cumulants of any given order. 

To prove Eq. (fT5|) . we compare Eq. (fT3|) with the defi- 
nition of the cumulant generating function 



io g % l (z)^j2—r K 

3 = 1 J ' 



(16) 



to deduce the remarkable identity 



ay(z) — n(n + v) + 



try 



(17) 



Substituting it back to Eq. (|M|) . we discover Eq. (|15|) as 
well as the above stated initial conditions. 

Large-n limit of the theory. — The nonpeturbative solu- 
tion Eq. (fTij|) has a drawback: it does not supply much de- 
sired explicit dependence of conductance cumulants Kj's 
on j. To probe the latter, we turn to the large-n limit of 
the recurrence Eq. (fT5|) . For simplicity, the asymmetry 
parameter v will be set to zero. 

Since, in the limit of a large number of propagating 
modes (n ^> 1), the conductance distribution is expected 
13711 to follow the Gaussian law 



i 



v /27rvar 00 (G) 



exp 



(.9 - n/2) 2 
2 varoo(G) 



(18) 



with the average conductance E[G] = n/2 and the con- 
ductance variance var^G) = 1/16, it is natural to 
seek a large-n solution to Eq. (|T5|) in the form kj — 
(rt/2)(5j,i + (l/16)(5 Ji 2 + Snj, where 5nj (with j > 3) 
account for deviations from the Gaussian distribution. 
Next, we put forward the large-n ansatz 



1 



E 



a m (j) 



n 1 L — ' n" 

m=0 



(19) 



which, after its substitution into the recurrence, yields 
the explicit formula 



_ 1 (2j - 1)! 



1 



m 2 i) 

8n 2 



(20) 



All odd order cumulants vanish identically. 



Interestingly, Eq. ([2H)) makes it possible to analytically 
study a deviation of conductance distribution f n {g) from 
the Gaussian law (g) . The Gram-Charlier expansion 



fn(g) = CXp 




(21) 



is the key. As soon as \d g log (g) \ ~ n, the operator in 
the exponent is dominated by the m = term in Eq. ((19)) . 
This observation reduces Eq. (|21[) to 



/«(<?) 



d\e~ n A 



2n^ 4 



r(i/8) V 7T y A 7 /8VTT2A 



exp 



2n 2 n 2 
1 + 2A 



(22) 

Here, r\ is the rescaled conductance r\ = 2(g/n) — 1. 

Equation (|22p is particularly suitable for the asymp- 
totic analysis. Performed with a logarithmic accuracy, it 
brings: 



log fn(g) 



-2n r] , 



\f]\ < 



1 



-2n 2 



l«|-^-^ lo g^, \ < \v\ < 1 

(23) 

This result shows that the Gaussian approximation 
for the conductance distribution is only valid for 
\g — n/2\ < n/4. Away from this region, the conduc- 
tance distribution exhibits long tails described by the 
exponential rather than the Gaussian law. Finally, it is 
straightforward to derive from the Toda lattice Eq. (fT0|) 
that, in the vicinity \g — g*\ < 1 of the edges [HI g* = 
and <?* = n, the conductance distribution exhibits even 
slower, power-law decay 2(| 23 1 



log/nfc) ~ (n 2 -l)log(2\ri-r],\)- — + — log n (24) 



with n* = ±1. 

Conclusions. — We have shown that a marriage be- 
tween the scattering S-matrix approach and the the- 
ory of integrable systems brings out an efficient for- 
malism tailor-made to analysis of the universal aspects 
of quantum transport in chaotic systems with broken 
time-reversal symmetry. Having chosen the paradigmatic 
problem of conductance fluctuations in chaotic cavities 
with ideal leads as an illustrative example, we determined 
the cumulants of conductance as well as its distribution 
exactly for any given number of propagating modes in the 
leads. It should be stressed that the ideas presented in 
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the Letter can equally be utilised [27[ to describe statis- 
tical properties of the shot noise power and the dynamics 
of charge transfer. 

Certainly, more effort is needed to accomplish inte- 
grable theory of the universal quantum transport. Ex- 
tension of the formalism presented to the (3=1 and 4 
symmetry classes and waiving the uniformity of the §- 
matrix distribution are the two most challenging prob- 
lems whose solution is very much called for. 

This work was supported by the Israel Science Foun- 
dation through the grant No 286/04. 

Note added. — Recently, we learnt about the paper by 
M. Novaes [38| who noticed that the n- th moment of 
conductance can nonperturbatively be calculated by us- 
ing the machinery of hypergeometric functions of matrix 
argument. Neither Toda lattice nor Painleve V represen- 
tations for the conductance distribution surfaced there. 
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